Effective temperature of active complex matter 
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We use molecular dynamics simulations to 
study the dynamics of an ensemble of interact- 
ing self-propelled semi-flexible polymers in con- 
tact with a thermal bath. Our intention is to 
model complex systems of biological interest. We 
find that an effective temperature allows one to 
rationalize the out of equilibrium dynamics of the 
system. This parameter is measured in several in- 
dependent ways from fluctuation-dissipation re- 
lations and by using tracer particles and they all 
yield equivalent results. The effective tempera- 
ture takes a higher value than the temperature of 
the bath when the effect of the motors is not cor- 
related with the structural rearrangements they 
induce. We show how to use this concept to in- 
terpret experimental results and suggest possible 
innovative research directions. 

Introduction.- Active matter entities, be them parti- 
cles, filaments or other, absorb energy from their envi- 
ronment or internal fuel tanks and use it to carry out 
motion. In this ubiquitous type of soft condensed mat- 
ter, energy is partially transformed into mechanical work 
and partially dissipated in the form of heat Units 
interact directly or through disturbances propagated in 
the medium. The effect of motors can be dictated by 
the state of the unit and/or its immediate neighborhood 
and it is not necessarily fixed by an external field. Active 
matter is thus kept in a non-equilibrium steady state and 
presents a number of interesting dynamical features, un- 
usual mechanical properties, very large responses to small 
perturbations, and very large fluctuations. Realizations 
of active matter in biology are manifold and exist at dif- 
ferent scales. 

Surprisingly enough, studies of active matter have been 
mainly analytical so far, based on refined calculations us- 
ing very stylized models 0-0, computer simulations of 
continuum models [l(J, and numerical studies of rela- 
tively simple lattice models [TJ E3l • Molecular dynamics 
simulations, that have proven to be so helpful to eluci- 
date the behavior of complex passive systems, have not 
been much employed in this field yet. 

Quite independently, passive systems with complex out 
of equilibrium dynamics have been the focus of much 
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study. These are systems with competing interactions 
between their constituents and no external energy input 
that in some conditions (such as sufficiently low temper- 
ature or high density) cannot equilibrate with their en- 
vironment. Glassy systems are the typical example. An 
important outcome of the theoretical modeling of these 
systems, later confirmed with numerical simulations, is 
the generation of an effective temperature, T e ff [1 31 ] - 
T e ff should be understood as a time-scales dependent 
parameter which takes a thermodynamic meaning in sys- 
tems with slow dynamics only [3, H3 |. Experimental 
measurements of T e f f in different glassy systems [l5l - [20j 
have been performed. 

In [2l| we took a first step in the direction of analyz- 
ing T e ff in active matter. We studied the dynamics of an 
ensemble of self-propelled particles, meant to be a rea- 
sonable model for real but simple active matter such as 
bacterial colonies, with molecular dynamics simulations. 
Here we get closer to more complex cases and we study 
a model of filamentous semi-flexible polymers [22|. We 
analyze its structure and dynamics and we study T e ff 
by using a variety of independent measurements. In par- 
ticular, we demonstrate that tracer particle techniques, 
that have already been exploited in the study of the me- 
chanical properties of real biological systems 0, [24j , can 
also be useful to obtain a direct characterization of the 
out-of-equilibrium state. Details of the calculations will 
be given in a forthcoming publication (25j . 

The model.- We consider the model for a passive semi- 
flexible linear polymers (filaments) of Ref. [26J . The opti- 
mal values of the parameters have been fixed by prelim- 
inary calculations [25|, [2(|. Polymer chains are coarse- 
grained and each segment is formed by identical beads, 
with mass m and diameter tq. Each bead is in contact 
with a thermal environment that is described by a ran- 
dom noise and a viscous drag. The dynamical evolu- 
tion at a bath temperature T is therefore controlled by a 
Langevin equation. 

The deterministic mechanical conservative force 
among beads is given by intramolecular and inter- 
molecular contributions. Connectivity of the chains is 
assured by harmonic springs acting between nearest- 
neighbor monomers. Chain rigidity is controlled by 
a purely repulsive interaction between next nearest- 
neighbor monomers. The actual choice of parameters 
allows one to control the flexibility of the chains, ranging 
from semi-rigid to completely flexible '2$, 26| . 

Having in mind that no simple technique can take into 
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FIG. 1: (Color online.) Probability distributions of the gyra- 
tion radius, R g , at the investigated motor activities, /. The 
structure is highly non trivial with an average value shifting 
towards lower values upon increasing /. The passive case 
(/ = 0) is also shown. 



account chemical activity in a molecular dynamics cal- 
culation, physical insight has been used to choose a rea- 
sonable model for motor action. Only a fixed fraction 
of polymers are provided with motors that have a gentle 
but still detectable effect on the behavior of the system. 
We localize the non-conservative motor activity on the 
monomer at the center of the polymer (27j . The motor 
effect is given by a time series of random isotropic kicks 
generated by a suitable stochastic process [2lj. The ac- 
tion of the motors is independent of the structural re- 
arrangements they induce (adamant [28] motors). The 
force of all motors, / , is the same and a fraction of the 
mean conservative mechanical force acting on the equiv- 
alent passive system, F. We thus quantify the motor 
activity with the control parameter / = f M IF. More 
details on the model are given in Refs. [2lL [25f| . 

We express all quantities in dimensionless reduced 
units [2|| that we associate to reasonable values of the 
experimentally measured counterparts [3(|. The typical 
unit of energy is 2/cbT; typical lengths are measured in 
units of 0.4 nm implying forces of order 20 pN at ambient 
temperature, T = 300 K. 

Structure.- We first investigate the structure of the 
semi-flexible active polymer melt by keeping all external 
macroscopic conditions fixed and varying the motor ac- 
tivity /. The analysis of the static structure factors and 
pair distribution functions (not shown) is consistent with 
a picture in which the average nearest neighbor distance 
between beads pertaining to different polymers decreases 
upon increasing /. Overall, the motors have the mixed 
effect of making the melt more compact and more dis- 
ordered simultaneously. This picture is corroborated by 
the extremely complex probability distributions of the 
gyration radius, R g , which are highly non-Gaussian and 
very asymmetric (Fig.[IJ. Motor activity has the effect of 
pushing the filaments closer and these, at the same time, 
fold substantially. We interpret this behavior in terms of 
a competition between two effects: crowding, which is re- 
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FIG. 2: (Color online.) a) Time-dependence of the inter- 
mediate scattering function, F 3 (Q,t), at the wave-vector of 
the first maximum of the static structure factor (activities 
increase to the left), b) Diffusion coefficient, D, and inverse 
structural (self-)relaxation time, r a , as a function of /. Data 
are normalized to the value for the passive system and follow 
power laws (dashed lines), with the indicated exponents. 



lated to excluded volume effects and tends to slow down 
the dynamics, and folding of the chains, which dynami- 
cally lifts the topological constraints, therefore decreas- 
ing entanglement. These have important consequences 
on the long-time evolution of the system, as it is verified 
below. 

Dynamics.- We next quantify the effect of motor ac- 
tivity on the long time dynamics. Figure [2] a) shows 
the self-intermediate scattering function, F S (Q, t), at the 
wave- vector corresponding to the main peak in the static 
structure factor and as a function of time for the inves- 
tigated motor forces. The mean-squared displacement 
data (not shown) confirm the picture whereby the collec- 
tive dynamics of the filaments gets faster under stronger 
motors. This is consistent with the gyration radius data 
and supports our interpretation: the augmented fold- 
ing of the filaments seems to solve local topological con- 
straints and diminishes the effect of entanglement. A 
quantitative analysis is presented in Fig. [5]b), where we 
plot the dependence of the diffusion coefficient and the 
relaxation time upon the strength of the motors. Quan- 
tities are shown renormalized to the value of the pas- 
sive system. It is very interesting to note that the two 
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FIG. 3: (Color online.) Fluctuation-Dissipation relations at 
the indicated values of /. The red data on top (closed circles) 
are for the passive system and the dashed (red) line has slope 
— 1/T confirming that FDT holds in this case. The solid lines 
are linear fits used to estimate the value of T e //. 



sets of data are consistent with very simple functional 
forms (dashed lines), D/D f=0 ~ 1 + 1423 x f 2 W and 
(t/t/ =0 ) _1 — 1 + 19 x /, respectively. 

Strikingly, the /-dependence of the diffusion coefficient 
is reminiscent of what was found in [3l], [32|. For the 
active colloidal particles used in these studies, the ef- 
fective diffusion coefficient can be written as D/D a = 
(1 + l/9P e 2 ), where D Q is the diffusion coefficient of the 
passive system and P e is the Peclet number, which char- 
acterizes the particle activity [U, |33j]. In our case the 
exponent (2.3 ± 0.1) is close to 2, the differences being 
possibly ascribed to the uncertainties on the determina- 
tion of D or to the presence of corrections to the proposed 
scaling. We speculate that here / plays the role of the 
Peclet number. A similar rationale should be valid for the 
striking linear dependence of the inverse of the relaxation 
time. 

Effective temperatures. - The analysis above proves that 
both structure and dynamics of an active system are in- 
fluenced by motor activity in a very complex fashion. 
We ask now whether it is possible to embed all this com- 
plexity in a single parameter, also prone to be directly 
determined in experiments. The chosen quantity is the 
effective temperature T e ff. One way of measuring the 
effective temperature in a non-equilibrium steady state 
consists in using the generalization of the fluctuation- 
dissipation theorem to out-of-equilibrium conditions [l3[ . 
The most visual way of analyzing the data is to con- 
struct a parametric relation between the integrated lin- 
ear response, XAB(t), and the correlation Cab{1) of a 
wisely chosen pair of observables, A and B, and asso- 
ciating minus the inverse of its slope with a possibly 
time-dependent parameter T e ff(t) [HI, H3j. Such mea- 
surements were performed in different out-of-equilibrium 
conditions [35l438l | demonstrating that T e f f (t) is (asymp- 
totically) piecewise constant taking, typically, two values: 
the one of the environmental bath at short times and an 



additional one characterizing the structural relaxation at 
long times, in cases with slow dynamics. 

We performed these measurements for different values 
of the forcing / and the results are shown in Fig.[3J (The 
needed observables are chosen such that the considered 
correlation functions are the F S (Q, t) of Figure [2] a). De- 
tails can be found in Ref. [Hj].) The curves clearly show 
two regimes. The first one, corresponding to high fre- 
quencies, represents the fast vibrations of the monomers, 
before any structural relaxation takes place. The second 
regime, at long times, describes the actual structural re- 
laxation in the sample. In all cases, the parametric plot 
is rather well described by a straight line, implying that 
T e ff is a constant over this time-scale. These data are 
also shown in Fig. 3] with open (red) squares. 
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FIG. 4: (Color online.) Summary of the entire set of data 
generated in this study. The power-law dependence of the 
effective temperature on motor activity, T e f / /T = 1 + jf 2 , is 
also shown (dashed line). 

The construction detailed above is very powerful, di- 
rect, and easily exploitable in computer simulations. In 
experiments, however, its implementation is far from ob- 
vious. This difficulty has been circumvented by consid- 
ering the dynamics of tracers (or intruders) , normally of 
spherical shape. These objects of micrometric size are 
nowadays quite simple to control and manipulate by ac- 
tive microrheology techniques. When immersed in an 
extended complex environment, the latter acts as an ex- 
ternal reservoir and the intruders' dynamics provide a 
wealth of information on the underlying matrix state. In 
our simulations we have implemented these ideas in two 
ways. 

A possible realization of an adapted thermometer that 
couples to the structural relaxation of the sample [l3[ is 
a free tracer particle with mass much larger than the one 
of the polymer beads. Its kinetic energy can be related 
to the system's effective temperature via an equipartition 
theorem (3(| [Hj|. Our results are shown in Fig.|4]as (blue) 
open circles. 

Another measurement of T e f f is based on the com- 
parison between (free) diffusion and (driven) mobility of 
tracer particles immersed in the active sample. Conform- 
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ing to the Einstein relation, T e f t should be the parameter 
linking the tracers' mean-square displacement and their 
driven displacement under a weak external applied force, 
in the long time limit in which the dynamics is diffusive. 
Data for T e ff obtained in this way are included in Fig. [4] 
as open (green) triangles. 

Figure @] collects all our results in the form T e ///T 
against /. We also include (diamonds) average values 
with error bars calculated from the entire set of data. The 
figure demonstrates that independent determinations of 
T e ff yield consistent results and suggests that, in our 
system, there is a single frequency-independent T e f f pa- 
rameter, with a precise thermodynamic meaning. Data 
tend to T e ff — T in the limit of vanishing activation, 
as they should, and the deviation from the environmen- 
tal temperature monotonically increases with increasing 
forcing. Most importantly, data are represented by the 
empirical law T e ff/T = 1 +7/ 2 , with 7 = 15.41 (dashed 
line). Remarkably, this finding supports our conjecture 
that the parameter / plays here a role analogous to the 
Peclet number for colloidal active particles considered in 
a few very recent experiments (3ll. |33|| . 



In conclusion, we studied the out of equilibrium be- 
havior of a system of active motorized filaments, a model 
for complex biological structures, using a genuine mi- 
croscopic approach that is still quite innovative in the 
field. We described the influence of motor activity on 
both structure and dynamics, finding a rationale for very 
similar recent experimental observations. We confirmed 
the relevance of the effective temperature, a concept well 
established in statistical mechanics of glassy systems, 
for complex active matter. Most important, we demon- 
strated that this parameter, relatively simple to measure 
by microrehology techniques, can be directly related to 
chemical activity, a process of extraordinary complexity. 
It would be interesting to explore the possibilities offered 
by the molecular dynamics approach used in this paper 
to clarify the extremely complex phenomenology of cell 
mechanical stability [HI, |4(j. In particular, one could try 
to give an answer to the question as to how mechanical 
properties (such as elastic moduli) change with motor ac- 
tivity and whether these changes can be rationalized in 
terms of an effective temperature. 
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